International Journal of Engineering Research and Development 

e-ISSN: 2278-067X, p-ISSN: 2278-800X, www.ijerd.com 
Volume 8, Issue 8 (September 2013), PP. 06-12 



Water Wave Generation Due to Initial Disturbance in 
Presence of an Inertial Surface in an Ocean with Porous 

Bed 

Subhabrata Gangopadhyay 1 , Uma Basu 2 

department of Mathematics, The Bhawanipur Eduction Society College, 5,L.L.R. Sarani, Kolkata- 700 009, 

India 

Department of Applied Mathematics, University of Calcutta, 92, A.P.C Road, Kolkata - 700 009, India 



Abstract:- The proposed research paper is concerned with the generation of water waves due to initial 
disturbances at the upper surface of a fluid assumed to be covered by an inertial surface. The inertial surface is 
composed of thin but uniform distribution of non -interacting floating material. The fluid is assumed to be of 
uniform finite depth but the bed is considered to be porous. The Fourier and Laplace transform techniques have 
been used as the mathematical analysis to obtain the depression of the inertial surface in the form of an infinite 
integral. The inertial surface depression is evaluated asymptotically for large time and distance using method of 
stationary phase for an initial disturbance concentrated at a point. They are also depicted graphically and 
appropriate conclusions are drawn. 
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I. INTRODUCTION 

The classical two-dimensional problem of water wave generation due to an initial disturbance in the 
free surface of the ocean was studied in treatises of Lamb[l] and Stoker [2] using linear theory of water waves. 
They used the method of Fourier transform and obtained the free surface elevation in the form of infinite 
integrals which were then evaluated approximately for large time and distance from the source of disturbance. 

This problem is referred to as Cauchy-Poisson problem in the water wave literature. 
Kranzer and Keller [3] discussed the problem of water wave generated by explosions in case of three 
dimensional motion in water of finite depth due to an initial surface impulse or initial surface elevation on a 
circular area. They also compared their theory with the available experimental results. Chaudhuri[4] extended 
these results in case of any initial surface impulse and elevation across arbitrary regions. In this respect, Wen[5] 
derived asymptotic results using the method of stationary phase as applied to double integrals. Debnath and 
Guha[6] worked out asymptotic analysis to the integral solutions while investigating the form of the axially 
symmetric free surface response of an inviscid stratified fluid in presence of an initial displacement of the free 
surface of the fluid of finite or infinite depth. Mandal[7] considered water wave generated by disturbance at an 
inertial surface. He employed Laplace transform technique to solve the initial value problem describing waves 
generated by a disturbance created at the surface of water covered by an inertial surface composed of a thin but 
uniform distribution of floating particles. Ghosh et.al[8] examined water wave due to initial disturbance at the 
inertial surface in a stratified fluid of finite depth. Dolai[9] considered generation of interface waves due to 
initial disturbances in the form of an initial depression or impulse at the interface when the upper fluid extends 
infinitely upwards and the lower fluid extends infinitely downwards. Dolai and Banerjea[10] considered the 
problem of generation of capillary-gravity waves due to initial interface disturbance at the interface between two 
superposed fluids, in which the upper fluid extends infinitely upwards and the lower fluid extends infinitely 
downwards. Ghosh[ll] considered generation of interface waves due to initial axially symmetric disturbances at 
the inertial surface between two superposed fluids wherein the lower fluid is of uniform finite depth while the 
upper fluid extends infinitely upwards. Maity and Mandal [12,13] threw considerable light on the problem of 
generation of surface waves caused by an initial disturbance at the upper surface of the ocean with an elastic ice- 
cover. They applied Laplace-Fourier transform to solve the problem and derived the elevation of the upper 
surface and also computed numerical results leading to graphs depicting the variations of surface elevation for 
large values of time and distance. In the present paper, the Cauchy-Poisson problem is considered, the upper 
surface of the ocean being covered by an inertial surface while the bed of the ocean is assumed to be porous in 
nature. 

Study of problems involving inertial surface has gained importance as sometimes in cold regions the 
ocean is covered by broken ice, which constitutes an inertial surface. 
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Fourier and Laplace transform techniques are employed to solve the problem. After invoking the inverse 
transforms, the potential function describing the wave motion is obtained and then the form of the inertial 
surface is deduced in terms of an infinite integral. This is then evaluated asymptotically for a large value of time 
and distance by using the method of stationary phase when the initial disturbance at the inertial surface is 
concentrated at the origin. The asymptotic form of the inertial surface is graphically presented. The figures 
depict variations of the non-dimensional inertial surface depression at a fixed point x f for different time t \ and at 
a fixed time t' for different x\ x f and t' being the non-dimensional horizontal distance from the origin and time 
respectively. 

Plots are considered for various values of the non-dimensionalized inertial surface parameter E f and for 
different values of the non-dimensionalized porous parameter G'h. 

The novelty involved in the current paper lies in the fact that the natural seepage of the ocean bed is 
taken into account as the bottom of the ocean is considered to be permeable unlike the previous authors who 
considered the bed to be impermeable. 

II. MATHEMATICAL FORMULATION OF THE PROBLEM 

We consider the upper surface of the fluid to be consisting of an inertial surface composed of non- 
interacting floating material of area density pE where p is the density of the liquid and E = corresponds to a 
liquid with a free surface. The ocean is assumed to be of uniform finite depth h but the bed is assumed to be 

porous where G ' = ^= is the porous effect parameter. The quantity a is a dimensionless constant which depends 

on the structure of the porous medium and v is permeability of the porous medium. 

The motion of the fluid is generated due to an initial disturbance prescribed on the inertial surface in 
the form of an initial depression of the inertial surface. We choose a rectangular co-ordinate system such that 
y = denotes the undisturbed inertial surface and y = h denotes the ocean bed, the y-axis pointing vertically 
downwards. Two-dimensional motion is considered and the motion of the fluid remains irrotational since it is 
assumed to start from rest. The velocity potential cp(x,y)t) describing the motion of the fluid satisfies the 
following initial value problem given by: 
V 2 cp = 0,0 < y < h,t > (1) 

[cp - Ecp y ) - gcp y = on y = 0, t > (2) 

j--G'<p = Oony = h,t > (3) 
where g is the acceleration due to gravity. 

The initial conditions when the initial disturbance is in the form of a prescribed depression of the inertial surface 
are: 

cp - Ecpy = at t = for y = (4) 

— {(p- E(p y ) = gf{x) at t = for y = (5) 

where /(x) is the prescribed initial depression of the inertial surface. 

III. METHOD OF SOLUTION 

Let cp(k, y; t) denote the Fourier transform of cp(x, y; t) given by cp(k, y; t) = -= J_ cp(x, y; t)e~ lkx dx 

On applying Fourier transform, equations (1) - (5) reduce to 
(p yy - k 2 (p = 0, 0<y<h f t>0 \ 



[cp — Ecpy) — gcpy = on y = 0, t > 
— G cp = ony = h,t > 



dt 



(6) 



dy 

cp — Ecpy = ony = at t = 

— {cp - Ecpy) = gf(k) on y = at t = 7 

where / (/c) denote the Fourier transform of f(x). 

Next, let cp(k, y; p) denote the Laplace transform in time defined by 

cp(k,y;p) = cp(k,y;t)e~ pt dt, p > 

Applying Laplace transform in time, equations in (6) reduce to 
y yy -k 2 y = 0,0 <y<h (7) 

v 2 y -(g + Ep 2 )v y = gf(k) on y = (8) 

^- Gy = on y = h (9) 
Equation (7) suggests: 
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B = A x e~ ky + A 2 e ky ; ^ = -A x ke~ ky + A 2 ke ky 

dy 

Using conditions (8) and (9) we obtain: 
(p 2 + kg)A 1 + A 2 p 2 (l - Ek) = gf{k) 

(ID 



(10) 



A 1 k-G' 

From (10) and (11) we solve for A 1 and A 2 and get: 

gfm i 



A, = 



B(k) p 2 +A 2 



(12) 



A 2 = k+G ' c -2kh . 9f(k) 



k-G' 

where A 2 = — 



B(k) p 2 +A 2 

(14) 



(13) 



and B(k) = 1 + (1 - Ek) ■ — e - 2kh 

k G 

Therefore we have, 



B(k) 



-ky 



p 2 +A : 



+ 



k+G _ r -2kh 
k-G' 



(15) 



,ky 



(16) 



If the initial displacement of the inertial surface is assumed to be concentrated at the origin, then we can take 

- \ 

f (x) = 8(x) or, / (/c) = -= where 8(x) is Dirac's delta function. 

y2n 



9 



Hence, (p(k,y;p) = - 

>vv >s>r; ^2nB(k) L k-G' 

Taking inverse Laplace transform we get, 



ky + e -2fc/i . e /cy 



p 2 +A' 



(17) 



(p(k,y; t) = 



9 



[e~ ky 



V2ttF(/c) L " k-G' 

The Fourier inversion then yields 



(18) 



. e ~2kh . gfey 



sinAtcoskxdk 



(19) 



The depression of the inertial surface at time t is then given by 
rj(x,t) = -<p t (x,0,t) 

9 



1 f 00 1 L , k+G' -ohh 



k-G' 



cosAt coskx dk 



= - J °° C(k)cosAt coskx dk, (20) 
where £(/c)C(/c) = \ + k -^e~ 2kh (21) 



k-G' 



IV. ASYMPTOTIC ANALYSIS : STATIONARY PHASE METHOD 

Equation (20) is an oscillatory integral which can be approximated by employing the method of 

X 

stationary phase for large x and t such that - remains finite. 
We can write, r)(x, t) = ^ Re[f™ C(k){e lt s ^ + e lt s ^}dk (22) 



where, S(k) = A — k ■ - = 



kx 
t 



(23) 



The second term of the integral (22) has no stationary point within the range of integration. The 
stationary point of the first integral is given by the following equation: 

(24) 

In view of the fact that the function S (/c) is monotone decreasing in nature, the stationary point k is 
the unique positive real root of the transcendental equation (24) and the integral (20) representing the free 
surface elevation can be approximated according to the stationary phase method as: 



r]{x,t) = C(fc ) 



]2nt \s\k Q )\ 



COS 



(25) 



V. GRAPHS 

We have considered the case when the initial disturbance is applied in the form of initial depression at 
the origin by taking f{x) = 8{x). In this case, the approximate form of the free surface depression is given by 
asymptotic expression (25). 

The physical quantities are non-dimensionalized as follows: 



x y t 

x =p y =v t = 
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To depict the effect of the inertial surface and the bottom porosity on the wave generation by the initial 
disturbance, the non-dimensional asymptotic forms of rj'(x'* are displayed graphically against x f for fixed t' 
and against t' for fixed x f in a number of figures. In order to perform it, the unique positive root of S (/c) = is 
obtained first. In view of the fact that the function S (/c) is monotone decreasing in nature, such a root exists 
uniquely. Clearly, k is a function of x' and t' and other parameters and these are computed numerically for 
different values of x ', t f and other parameters. 

To visualize the nature of the wave motion generated, r\'(x , t ) are plotted in Figs.l, 2 against x' with 
fixed t ' = 20 where in Fig.l, G h = and in Fig. 2, G h = 0.1; 77 (x , are plotted in Figs. 3, 4 against x' with 
fixed t = 30 where in Fig. 3, G h = and in Fig.4, G h = 0.1; in each of the Figs. 1,2,3,4 E = 0, 0.1, 0.01 and 
p = 0.2. From Figs. 1,2,3,4 it is observed that for fixed t\ when x f increases, amplitude of the wave profiles 
decreases so that at large distances they die out. This seems feasible as the initial disturbance is concentrated at 
the origin. From Figs. 1,2,3,4 it is also observed that in absence of porosity i.e., when G h = 0, the curves 
corresponding to E' = 0, 0.1, 0.01 almost coincide while at G h = 0.1 the effect of differences in the values of 
the inertial surface parameter is more pronounced. This is possibly due to an interaction between porous 
parameter effect arriving upwards from the bottom and the inertial surface parameter effect arriving downwards 
from the surface of the water. 

In Figs. 5,6 rj'(x ,t) are plotted against t f with fixed x ' = 10 where in Fig. 5, G h = and in Fig. 6, 
G h = 0.1; r]\x , t ) are plotted in Figs. 7,8 against t' with fixed x = 15 where in Fig. 7, G h = and in Fig. 8, 
G h = 0.1; in each of the Figs. 5,6,7,8 E = 0, 0.1, 0.01 and p = 0.2. From Figs. 5,6,7,8 it is observed that for 
fixed x\ when t f increases, amplitude of the wave motion at the inertial surface increases. This is somewhat 
unrealistic and is a consequence of the presence of a strong singularity at the origin as has been remarked upon 
by Stoker[2]. It is further observed from Figs. 5,6,7,8 that there have been a steady increase in the amplitude of 
the wave profile for E = 0.01, 0.1, respectively in all cases. 

All the eight figures demonstrate the interaction of the porous bed effect with the floating materials on the 
surface and both the phenomena affect the wave motion significantly. 




Fig.l Wave profiles of an inertial surface depression due to an initial disturbance in the form of an initial 

depression 

" ■ {=20.601=0.1 




Fig.2 Wave profiles of an inertial surface depression due to an initial disturbance in the form of an initial 

depression 
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Fig.3 Wave profiles of an inertial surface depression due to an initial disturbance in the form of an initial 

depression 
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Fig.4 Wave profiles of an inertial surface depression due to an initial disturbance in the form of an initial 

depression 




in 



Fig.5 Wave profiles of an inertial surface depression due to an initial disturbance in the form of an initial 

depression 
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Fig.6 Wave profiles of an inertial surface depression due to an initial disturbance in the form of an initial 

depression 
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Fig.7 Wave profiles of an inertial surface depression due to an initial disturbance in the form of an initial 

depression 
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Fig.8 Wave profiles of an inertial surface depression due to an initial disturbance in the form of an initial 

depression 

VI. CONCLUSIONS 

The problem of water wave generation due to an initial disturbance in the inertial surface of the ocean 
with uniform finite depth is investigated using linear theory of water waves. Fourier and Laplace transform 
techniques are used to solve the problem of the initial disturbance with an inertial surface. The asymptotic form 
of the non-dimensionalised inertial surface depression is evaluated for a large value of non-dimensionalised time 
and distance by using the method of stationary phase. Plots of rj (x ,t ) in various cases indicate the possible 
interaction between bottom porosity and the effect of floating material which composes the inertial surface. 
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